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(<-^ ' Abstract 

^^ ' We give a complete and rigorous proof of the Unruh effect, in the 

Oj. following form. We show that the state of a two-level system, uniformly 

I ' accelerated with proper acceleration a, and coupled to a scalar bose field 

i^H I initially in the Minkowski vacuum state will converge, asymptotically in 

j^ , the detector's proper time, to the Gibbs state at inverse temperature 

^ _ 22L The result also holds if the field and detector are initially in 
an excited state. We treat the problem as one of return to equilibrium, 
exploiting in particular that the Minkowski vacuum is a KMS state with 
L^ " respect to Lorentz boosts. We then use the recently developed spectral 

rS , techniques to prove the stated result. 

. P., 

1 Introduction 

The following observation, now referred to as the Unruh effect, was made by W. 
Unruh in 1976 23- When a detector, coupled to a relativistic quantum field in 
its vacuum state, is uniformly accelerated through Minkowski spacetime, with 
proper acceleration a, it registers a thermal black body radiation at temperature 
T = 2^^- This is the so-called Unruh temperature. In more anthropomor- 
phic terms jUWj . "for a free quantum field in its vacuum state in Minkowski 
spacetime M an observer with uniform acceleration a will feel that he is bathed 
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by a thermal distribution of quanta of the field at temperature T." This result 
has attracted a fair amount of attention, and generated considerable surprise 
and even some scepticism. For a review of various aspects of the subject, nice 
physical discussions of the phenomenon, and further references, we refer to |Ta| . 
|Wa| |FiiU| . The reason for the surprise is that, if you think of the vacuum 
as "empty space", then you will find it puzzling that a detector, accelerated or 
not, which may itself initially be in its ground state, will "see particles", since, 
after all, in the vacuum, there aren't any. In order not to be surprised, one has 
to remember that, of course, the vacuum is not "empty space", but the ground 
state of the field, and one should expect the detector to react to the presence of 
the field when it is accelerated through space. 

For example, if you were to drag a detector along a non-relativistic chain of 
oscillators in its ground state, you would certainly expect the coupling between 
the detector and the oscillators to excite both. The energy for this process is, in 
final analysis, furnished by the agent that drags the detector along the oscillator 
chain. 

What is nevertheless still surprising in connection with the Unruh effect is 
the claim that the detector "perceives" a thermal distribution of radiation at 
some particular temperature that only depends on the acceleration. To see 
what is precisely meant by these statements, it is helpful to get rid of the 
anthropomorphic terminology used above and in much of the literature as well 
as of all reference to particles or quanta, which turn out to be irrelevant to 
the discussion. This is what we will do below. It is worth pointing out in this 
connection that already in ^ , "detection of a particle" is defined by "excitation 
of the detector", and does therefore not presuppose the actual definition of what 
a particle precisely is, which is a tricky thing to do, as is well known (Fu) . In fact, 
the computations in the physics literature of the excitation probability of the 
detector can be seen to be perturbative computations of the asymptotic state of 
the detector (see |UW| for example). We therefore adopt the following simple 
formulation of the Unruh effect. Consider the coupled detector-field system. 
Suppose that initially it is in a product state with the field in the vacuum state. 
Now let the coupled system evolve. At some later (detector proper) time, the 
state of the system will no longer be a product state. Now trace out the field 
variables, to obtain the reduced state of the detector (which will be a mixed 
state, even if the initial state was pure). The Unruh effect states that the latter 
converges, asymptotically in the observer's proper time, to the Gibbs state at 
the aforementioned temperature T. Note that this is not by any means obvious: 
after all, a priori, it is not clear why the detector state should, asymptotically 
in time, converge at all, and even if it does, it is not obvious it should tend to a 
positive temperature state: a priori, it could have been any other mixed state. 

It is our goal in this paper to give a complete and rigorous proof of the above 
statement. The way we have formulated it makes it clear already that we think of 
it as a problem in the theory of open quantum systems in which a small system, 
here the detector, is coupled to a reservoir, here the field. Let us formulate our 
result somewhat more precisely. For a completely rigorous statement, we refer 
to Section 121 The model we consider is the one proposed in |UW| . which is 



itself a simplification of the model considered in pj] . The detector is modeled 
by a two-level system and the field is taken to be a massive or massless Klein- 
Gordon field. The observable algebra of the detector is therefore generated 
by "fermionic" creation/annihilation operators A, A^{a). The free Heisenberg 
evolution of the detector is A{a) = —iEA{a), where a is the detector's proper 
time. In other words, the free detector Hamiltonian is 

Hb = EA^A. 

The coupling between the field and the detector is realized via a monopole, and 
is ultraviolet regularized; it is sometimes referred to as a de Witt monopole 
detector (see [35] ). Suppose initially the detector-field system is in a product 
state ll)q with the detector in a state described by some density matrix p and the 
field in the Minkowski vacuum state. Let B be a detector observable and a'^{B) 
its Heisenberg evolution under the coupled dynamics, with coupling constant A. 
Then we prove that 

lim iOoia^AB)) = ^Tr e'^^^B + 0{X'). (1) 

Here /? = (kBT)-^ with T the Unruh temperature and Z/3,d = Tre-'^^^. The 
approach to the equilibrium state is exponentially fast. 

Our proof of this result is based on techniques developed in the last decade to 
prove "return to equilibrium" in open quantum systems |.TP1I I.TP2I IBaFSI IMll 
IDJ, DJPj . We combine these with the Bisognano- Wichman theorem |BiWi| , 
which states that the vacuum is a KMS state for the Lorentz boosts on the 
Rindler wedge. The relevance of this last result to the Unruh effect (and a 
generalization to more general spacetimes) was explained a long time ago by 
Sewell in [Sij. Let us point out that the work of Sewell, together with known 
stability results of KMS states (see e.g. jDal lKFGV| ) imply a result somewhat 
similar to but considerably weaker than Q, namely 

liin Ma'AB)) = ^Tr e-^"-B. (2) 

This is the so-called van Hove weak coupling limit. In our result, the limit 
cr ^ oo is shown to exist for all sufficiently small A, and to coincide with the 
right hand side of (QJ. 

The paper is organized as follows. In Section |2l we describe the model in 
detail and state our main result. We will also comment on the precise role 
played by the choice of the form factor determining the ultraviolet cutoff in 
the interaction term. Sectional is devoted to its proof. The latter uses Araki's 
perturbation theory for KMS states and its recent extensions, together with the 
spectral approach to the problem of return to equilibrium developed in the cited 
references. Since this material is rather technical, we have made an effort to 
state the result in Section 12 with as little reference to it as possible. 
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2 The model and the result 

We need to give a precise description of the model and in particular of its 
dynamics. This requires some preliminaries. 

2.1 The free field 

Let us start by describing in detail the field to which the detector will be cou- 
pled. The field operators are represented on the symmetric Fock space T over 
L^(R'^,dx). Here d > 1 is the dimension of space and x = {x'~',x) be is a 
point in Minkowski space-time M x M'' (with metric signature (+, — , . • • , — )). 
So T :— ®„eN-?^''"\ where JF'"-' is the n-fold symmetric tensor product of the 
one-particle space L'^ {M.'^ , dx) . 

Let 5(R''+^R) and 5(R''+^C) denote the real and the complex valued 
Schwartz functions on IR''+^, respectively. For / G 5(R''+^;C), one defines 
the field operators in the usual way: 

f] = (-A + m2)i/2, ^±/ = ^ dt-^e±'^*/t, Q[f] = -^(^HS+f) + a(S^)). 

Here A is the Laplacian, m > the mass, and a, a^ are the usual creation and 
annihilation operators on J^ (we follow the convention that / i-^ a^(/) is linear 
while / I— > a{f) is antilinear), and the bar denotes complex conjugation. When 
m = 0, we will suppose d > 1. Writing formally 

Q[/]=/ dxf{x)Qix), x^{x°,x), 

this leads to the familiar 

Q{x)= [ '^- pi£-i"(fc)^°»(^) + p-ifci£+i"(fc)^°„*m] (3) 



where w(fc) = y fc^ + m? . The field satsifies the Klein-Gordon equation \3Q{x) + 
w?Q{x) — 0, where D — d^o — A. We use units in which h — 1 — c. 

As can be learned in any book on special relativity (such as plj). in an 
adapted choice of inertial coordinate frame, a uniformly accelerated worldline 
of proper acceleration a > 0, parametrized by its proper time a, has the form 

X {a) = — sinhacT, x (a) = — cosh acr, x (cr) = = x (a), 
a a 



Associated to this worldline is the right wedge (or Rindler wedge) Wr := {x G 
R"*! x^ > la;"!}. It is the intersection of the causal future and past of the world- 
Hne, or the collection of spacetime points to which the observer on the worldHne 
can send signals and from which he can also receive signals. Note for later 
reference that the left wedge Wl := — Wr is the causal complement of Wu.- 

There exists a global coordinate system on Wb. that is particularly well 
adapted to the description of the problem at hand. It is given by the so-called 
Rindler coordinates {t,u,Xj_) eRxR*^ x R^~^, defined by 



usinhr, 



wcoshr, x^ = {x2, ■ ■ ■ Xd). 



(4) 



Here r is a global time coordinate on the right wedge. Note that, given a € 
R"*", (a2, • ■ • did) G M''"^, the curve u — 1/a, Xj_ = (a2, • ■ • , ctd) is the worldline 
of a uniformly accelerated observer with proper acceleration a and proper time 
a = a~^T. In addition, two points in the right wedge with the same value for 
the T-coordinate are considered as simultaneous in the instantaneous rest frame 
of any such observer (see ^ij). Among the Lorentz boosts, only the boosts in 
the xi -direction leave the right wedge invariant. In inertial coordinates they are 
given by the linear transformations 
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In the Rindler coordinates, this becomes Br'{T,u,Xj_) = {t + t' ,u,Xj_). In this 
sense, the boosts in the xi-direction act as time translations on the Rindler 
wedge. 

Since the field satisfies the Klein-Gordon equation, one has, in Rindler co- 
ordinates on Wr: 



{u-'^dl - u-^duudu + (-A_L + m^)) Q{t, u,x^) = 0. 



(5) 



Moreover, the covariance of the free field under the Poincare group yields, for 

all T e R, 

Q[/oi?_,]^c'^-^g[/]c-^-^ (6) 

with 

LY^dT{K), K^n^/'^X^^^/^, (7) 

where X^ is the operator of multiplication by x^ . In particular, for x = 
(r, u,x^) € H/R,, 

Q(r,u,x^) =e'^^"Q(0,u,a;^)e-'^^". 

In other words, Lp generates the free Heisenberg dynamics of the field operators 
associated to the right wedge. Let us furthermore introduce, for later purposes, 
the conjugate field 



P\I] ■■= -^Q[f ° B^r] \r=0 = i [ip, Q[/]] . 



(8) 



It then follows from the basic properties of the free field that the equal time 
commutation relations of the field and the conjugate field are, at t = 0, 

[QiO,u,x^),P{0,u',x:^)] = iuduiu') 4^ (ail). (9) 

The following useful identity follows from ij^l and ©: 

i[LF,P(0,u,x^)] = - {-uduudu + u^{-Aj_ + m^)) Q{0,u,xj_). (10) 

For an algebraic formulation of the dynamics, indispensable in what follows, 
we need to identify the observable algebra of the theory. The observable algebra 
of the field is ^f := {W{f)\f G 5(M''+\R)}", with W{f) = e-^^I^l the usual 
Weyl operators. One should think of the observable algebra as containing all 
bounded functions of the (smeared) field operators Q[f] or, more pictorially, all 
observables that can be constructed from the Q{x), x G R'*+^. Associated to 
the right and left wedges are local algebras of observables -4f;R,l := {T^(/)l/ G 
'5(Wr_l,]R)}". Again, those should be thought of as containing all observables 
that can be constructed with the field operators Q{x), for x belonging to the 
wedge considered. As pointed out above, one can define on Af an automorphism 
group a° by 

al^,{A) = e'^^^A e"'^^", A e ^f. (H) 

We note that a^ ^ leaves -4f;R invariant. 

2.2 The free detector 

As pointed out in the introduction, we think of the detector as a two-level 
system. Our results extend without problem to an TV-level system, at the cost 
of irrelevant notational complications. So we follow the physics literature on 
the subject and limit ourselves to a highly idealized two-level detector. Its 
observable algebra is simply the algebra of two by two matrices S(C^). It will be 
convenient to use a representation of this algebra in which both the ground state 
and the Gibbs state at inverse temperature (3 are represented by vectors. This 
representation, well known in the mathematical physics literature on quantum 
statistical mechanics, is of course different from the usual one in the standard 
physics literature in which the latter is represented by a density matrix. 

It is defined as follows. One represents the observable algebra S(C^) as 

Ad ■■= B{C^) (g) I2 on Hd = C^ (g) C^, with in particular ^t := q n ^ I2. 

The algebra ^d is generated by the identity operator, A\ A and A^A and one 
has AA"! + ylU = 1. 

In this representation, the free Heisenberg evolution of the detector with 
respect to its proper time a is generated by the self-adjoint operator Ld := 
-ffD® ]l2 - ]l2«'i?D, with 



Hd = 



E 




(12) 



for some E > 0, where E represents the excitation energy of the detector; Ld 
is referred to as the free Liouvillean of the detector. To see this it is enough to 
remark that 

satisfies the correct Heisenberg equation of motion 

A\a) = iEAHa) (13) 

of an unperturbed two-level system. Note that the energy levels of the detector 
are thought of in this model as pertaining to internal degrees of freedom f [Ul 
lUWLITaj l. One should think of the two- level system as being dragged through 
spacetime by an external agent that ensures it has constant acceleration a. So 
the translational degrees of freedom of the detector are not dynamical variables 
in this kind of model. In the representation above, the ground state of the 

detector can be represented by the vector | ) and the Gibbs state at inverse 

temperature P by the vector 

1/3, D) := (1 + c-^^)-i/2(| - -) + c-^^^/2| ^ ^^) ^ -j^^ 
Indeed, one easily checks that, for any B e 



(/3,D|S (g, l2|/3,D) = -^Tr e'^^^S, Z^.d = Tr e"''^°. 

It is the fact that both the ground state and positive temperature states of the 
detector can be represented by vectors that makes this representation particu- 
larly suitable for the problem at hand. 

2.3 The uncoupled field-detector system 

It is now easy to describe the observable algebra of the joint detector-field sys- 
tem, as well as its uncoupled dynamics. On the Hilbert space H :— Hb JF we 
consider the observable algebra A := Au ^ Af.r and the self-adjoint operator 
Lq — Ld (>^ly^ + t-Hn ^ aLp. The latter determines an automorphism group 

of A in the usual way: a°(B) = e'^o'^Be"'^"'^, B G A. Setting B{a) := a°{B) 
this yields a solution of the Heisenberg equations of motion of the uncoupled 
detector-field system on the Rindler wedge Wr, which are given by Q and il^ . 
with T = aa. 

We will be mostly interested in the state of the system where, initially, the 
detector is in its ground state, and the field in its Minkowski vacuum. This 

state is represented by the vector |g) := | ) |0) G H. We will write, for any 

Be A: 

(B), := (g|i?|g) (14) 



2.4 The coupled field-detector system 

For the coupled system we will use the same representation of the observable 
algebra, but change the dynamics. We will give a precise and mathematically 
rigorous definition of the dynamics below but to link it with the physics literature 
on the subject, we start with a formal computation. Let C{a) — [A{a),A'^{a)]. 
According to jUWj . the Heisenberg equations of motion of the observables of 
the coupled system are 

{U + m^)Q{x) = -\p{x,){A + A^)(^) (15) 

A{a) = —\EA{a)+i\C{(j) I dudx^au p{x*)Q{aa,u,XjJ, 

The function p tunes the coupling between the detector and the field. It is 
evaluated at 

x* := a; — x(r(x)/a), 

the spacelike vector linking x in the right wedge to the instantaneous position 
x{a) of the detector, at proper time a = T{x)/a, where t{x) is the Rindler time 
coordinate defined in Q. 

Let (r, u, Xj_) be the Rindler coordinates of the point x then the ones of x{a) 
are (r, l/a,0^) and hence we may identify a:*, whose coordinates are (0,u — 
l/a,Xj_), with an element of (— 1/a, +00) x M''^^. We take the coupling function 
p to be in C(f ((— 1/a, +cx)) x R''"^), normalized as / pix^dx"^ = 1. Typically we 
imagine p to be peaked at the origin, so that the field is coupled strongest at the 
position of the detector. Only for such couplings does it make sense to interpret 
a as the proper time of the detector. Indeed, if the detector is coupled to the 
field over a large spatial region, different parts of the detector undergo a different 
acceleration and have a different proper time. The mathematical result we 
obtain then still holds, but does no longer have the same physical interpretation. 
A coupling strictly localized at the position of the detector is formally given by 
p(x) — S{x), a situation which does not fit the rigorous mathematical setup 
presented in this work. We will comment further on the role played by the 
choice of coupling in Section [231 

Using and ljl()|l . it is easy to show through a formal computation that 
these equations are satisfied by the operators Q^^'{t, u,Xj_) and A'-^'^^a) defined 
as follows: 

Q^^H'r,u,x^) :=e'^^» Q{0,u,x^) e"'^^^, ^(^^(ct) := e'^^'^ Ae-'^''", 

where 

Lx := Lq + \I, I := {A + A'') / dudx^au p{x^\r=o)QiO, u, x^) (16) 



and Xi,\r=Q is given in Rindler coordinates by (0,u— l/a,a;_L). In other words, 
the Liouvillean L\ generates the correct Heisenberg dynamics of the observables 
in the representation at hand. 



Remark. The analysis we carry out in this paper works for general interac- 
tions of the form I — G ■ Q{g) + G* • Qi'g), and for sums of such terms, where G 
are matrices acting on the detector space, and g € L^(M^, dx) are "form factors". 

The following result is proved in Sectional 

Proposition 2.1 The operator L\ in \16\) is for all A essentially self-adjoint 
on D{Lo) n D{I) and the maps a^{B) := c^L^a^^-iL^a ^^^^ a & R and B & 
A define a weakly continuous one-parameter group of automorphisms of the 
observable algebra A. 

2.5 The result 

We are now in a position to give a precise statement of our result. Define 

g{>i, ki_) = {xip{x^ |r=o)) ((|fc± P + m2)i/2 sinh k, k^_) , (17) 

where ^ denotes the Fourier transform. 

Theorem 2.2 Let d > 1 if m > and d > 2 if m = 0, and suppose the 
following "Fermi Golden Rule Condition" holds, 

I dn e-'» ^ g[xAi_) + for some k^ G W^-^ . (18) 

Then there is a constant Aq > s.t. if < |A| < Aq then 

hm (a^(B))g = -^Tr e-^«°i? + 0{X'), (19) 

for all B e B{C^), and where (3 = ^. 

More generally, if g is any density matrix on H then 

lim Trga^{BF) = (- Tr e^'^^^fi) (0|i^|0) + 0{X^), (20) 

for any detector observable B e S(C^) and any field observable F G ^p- 

Result (I19II shows that if at ct = the detector-field system is in a state which 
is a local perturbation of its ground state, then the reduced density matrix of 
the detector converges asymptotically in time to the detector's Gibbs state at 
inverse temperature /3 = ^. This is a (slightly) stronger statement than the 
formulations usually found in the literature, since it allows both the field and 
the detector to be initially in an excited state. 

Remarks. 1) Theorem 12.21 follows from a more complete result, stated as 
Theorem l3.5l below. where the r.h.s. of Ill9|l is identified as the equilibrium state 
|A) S Ti of the coupled system, see also lf28jl below. An expansion of (A| • |A) for 
small A yields the uncoupled equihbrium state plus an error of second order in 



A (the absence of a first order error term is due to the fact that the expectation 
of the interaction / in the uncoupled equiHbrium state vanishes) . 
2) The approach to the Hmit state in ((Tflll is exponentially fast, 

\Trga^^iB)-{X\B\X)\<C\\B\\c-^"'^^, 

where C is a constant (depending on the interaction, but not on the initial 
density matrix g nor on B) and 77 = (1 + e~^'^^'°')£^, with 



^-m = ii dfc. 



"gi^iki 



> 0. (21) 



The quantity Tj-eiax = l/A^ry is called the relaxation time of the process. The 
purpose of condition (tT8|l is to ensure that ^ > 0, i.e., that Trciax < 00. We will 
show in the following subsection that this is typically the case. 

We finally remark that, whereas the leading term of the right hand side of 
iflfljl does not depend on the choice of form factor p in the interaction term, the 
relaxation time Trciax does, via ifTTjl and ll2Tll . Nevertheless, we show in the next 
subsection that Trciax is independent of the form factor for interactions sharply 
localized at the position of the detector. 

2.6 The Fermi Golden Rule Condition 

The goal of this section is to show that IJ18II is satisfied for "generic" interactions. 

Proposition 2.3 Take the coupling function p in U5\) . il6\) to be of the form 
pin) — Pi{xi)p±{x^), ("square detector") with pi > 0. Then condition ^^ is 
satisfied for all E except for E ^ E, where E is a discrete (possibly empty) subset 
ofR. In particular. S,{E) > for all E ^E. 

Values of E satisfying i,(Ei) — (which form necessarily a subset of E in the 
proposition) correspond to energy gaps of the detector Hamiltonian for which 
thermalization of the detector occurs (if at all) with a larger relaxation time at 
least of the order A"'' (as opposed to A~^ for E s.t. ^(E) > 0), see |M2| . 

For a particular choice of the coupling function p one may resort to a numer- 
ical study of the condition ltT8|l . On the analytic side we can calculate £,, (J2lJ, 
in the limit of a strictly localized interaction. More precisely, we choose pi, p± 
as in Proposition 1221 and consider the family Pe{x) = e~'^Pi{xi/£)p±{x^/e) -^ 
S{xi — l/a)(5(xj^) which represents an interaction localized exactly at the posi- 
tion of the detector in the limit e — > 0. Each e defines thus a Ce(-E) by lf2T|l . and 
we obtain, for d = 3 and tti > 0, 



Xxv^UE) 

e— tO 



a f dkj_ 
2 ~ 



2 sinh m: — cosh >c — \ 

t£i 

(-^-t-coshx)^ 



J LJ_L —11 — >fH — — sinh 3 
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where ijJ± — y/\kj_\'^ + m'^. This Hmit does not depend on the form of p and the 
leading term of Tj-ciax (as e ^ 0) is thus independent of the detector form factor. 

Proof of Pronosif,ion \2..'^ We denote the integral in ifTSJI by ip± (fc j^ )J{E,u!±), 
where uj± — y/\kj_\'^ + m'^, see also ifTTIl . For uj± ^ we can make the change 
of variable y = u!± sinh k to obtain the representation 

J{E,oj^) = / dy f{y), f{y) := e-'^/'^(-pi(y) + pi'(y)) . 

(22) 
We view lo\= [im the integral as a parameter, fi> 0. We first show that given 
any hq > 0, the integral in II22II . for E ^ 0, does not vanish identically in any 
neighbourhood of /io. 

Let us consider /io = 1; a simple modification of the following argument 
yields the general case. Assume ad absurdum that J(0, ^) = for all fj, in a 
neighbourhood of 1. Then, by taking derivaties of J{0,p) with respect to fJ,, at 
/Uq = 1, we see that 

"dy(l + y2)-"(l + y2)-i/2/(y)^0, (23) 

for all n = 0, 1, . . . Now, it is not difficult to verify that the linear span of ah 
functions (1 + j/^)^", n — 1,2, .. . is dense in the space of even functions in 
L^(M, dy). (One may prove this with little effort via the Fourier transform, 
for example.) It thus follows from i2'M that the even part of / must vanish, 
f{y) + f{—y) = for all y > 0. In particular, /(O) = 0, which means that 

a-i = -ipi'(O) (24) 

(we assume without loss of generality that pi is normafized as J^dx pi{x) — 1). 
On the other hand, we have — ipi'(O) — — J^dxxpi{x) < a~^, since in the 
integral, x > —a~^ due to the fact that pi is supported in (— l/a,oo). Therefore 
condition Ij24|l is not verified. 

This shows that given any /io > we can find a pi > (arbitrarily close to 
po) with the property that J(0,pi) ^ 0. 

Pick a nonzero Kq e M''"^ satisfying pi (/-Co) ^ and set po := ^/liiTop + m^. 
Then, by the above argument and by the continuity of pX there is a pi —: 
^\Ki\^ + mP- (which is close to po and defines a Ki close to Kq) s.t. J(0, pi) ^ 
andp±{Ki) ^ 0. Hence we have shown that there exists a nonzero K\ satisfying 
ipl(-ftri)J(0, wi) 7^ 0, where lo\ = •\/|^KiF+r7?. Condition ifTsll is thus satisfied 
for E = Q. 

Finally we pass to the other values of E by an analyticity argument. Indeed, 
one easily sees (best by using the form of J in which one integrates over >c 
rather than y, c.f. Ijl7|l . Ijl8|l l that the map E i-^ J{E,uJi) is analytic and by 
the previous argument it does not vanish at E = 0. Thus the zeroes of this map 
are contained in a discrete set £{uji) C C. Any E avoiding this set thus satisfies 
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3 Proof of Theorem EH 

3.1 Strategy 

As mentioned in the introduction, the first ingredient of the proof is the obser- 
vation that the Minkowski vacuum is (a reaHzation of) the GNS representation 
of a KMS state on the right wedge algebra for the Lorentz boosts at the inverse 
temperature f3 = 2tt. This is the content of Theorem 13.11 below. To give a 
precise statement, we need the so-called modular conjugation operator, defined 
as follows: 

-^F = r(jp), where V i/^ S L^{R'',dx), j^tp{x) = i^{~xi,Xj^); (25) 

here r(j) stands for the second quantization of j. 

Theorem 3.1 ([BiWiJ) The Fock vacuum in T induces on Af,r a state which 
is KMS at inverse temperature (3 = 27: for ap ^ . In particular, one has 



Apji—Af.l, JyAy.rJ^ — Ay,l, -4F,i?|0,F) — J^ — ylF,L|0,F), 
andfor all f eS{WK,C): iF|0,F) =0 and q-^ ^^ Q[f]\Q ,Y) = JfQ[7]|0, F). 

This result was proven in considerable generality in |BiWi| , for relativistic fields 
satisfying the Wightman axioms. The result above for the free scalar field can 
be obtained from essentially direct computations, and we shall not detail it. 

Similarly, the states |/?, D) introduced in Section [T2l are GNS representatives 
of the KMS states at inverse temperature (i for the free detector dynamics ao.cr 
on the detector observable algebra B{C'^). This well known observation is for 
convenience summarized in the following lemma. The appropriate conjugate 
operator is given by 

where C is the antilinear operator of complex conjugation on C^ and E is the 
exchange operator on C^ (Xi C^, Eip ®x — X®'-P- 

Lemma 3.2 For any /3 > 0, the vector \j3, D) induces on S(C^) a state that is 
KMS at inverse temperature (3 for a^^. In particular, one has 

A'j, = h^B{C^), A'd^JuAuJb, AB\P,D)=n = A'j,\P,D), 
and 

g-/3LD/2(^ ® 12)1/3, D) = Ju{B* (E) ]l2)|/3, D). 

Defining, on 7i = Hd ® ^, J := J^y ® Jy , it follows that the vector 

|0) :=!/? = 27r/a)(»|0,F) (26) 

is a GNS representative of the KMS state at inverse temperature (3 = — for 
the free dynamics a^^ on A = -4d ® Ay,r- This suggests to treat the problem 
at hand as one of return to equilibrium. 



12 



The rest of the argument then proceeds in three steps: 

(a) One proves the existence of a GNS representative |A) S 7i, defined below, 
of the KMS state for the perturbed dynamics at the same temperature (Section 
HI; 

(b) One reduces the proof of Theorem 13.51 and hence of Theorem 12.21 to 
showing that the generator of the perturbed dynamics has a simple eigenvalue 
at and otherwise absolutely continuous spectrum only; 

(c) One finally uses spectral deformation theory to prove these two state- 
ments. 

The strategy in (a)-(b)-(c) has been appHed successfully to radiative prob- 
lems in atomic physics, the spin-boson model, and similar systems in |JPlll,TP2l 
IBaFSI ImTI IDJPJ ■ where we refer for further references. A concise introduction 
to the field can be found in jPi]. The implementation of this strategy in the 
present context is reasonably straightforward. We will detail those points that 
are specific to the current situation. 

3.2 Perturbation theory 

We define on H = Hb ® ^, in addition to La (see El, the so-called standard 
Liouvillean ^ 

Lx = Lx- XJIJ. (27) 

We outline the proof of the following result in Section 

Lemma 3.3 L\ is essentially self-adjoint on D[Lq) n D{T) n D{JIJ) and, for 
all B eA, 

A useful feature of the standard Liouvillean (in fact, the motivation for its 
definition!) is that the unitary it generates leaves the equihbrium state of the 
coupled system invariant, see ll29ll below. 



Proposition 3.4 The vector |0) representing the uncoupled equilibrium state, 
\2f^) . is in the domain of the unbounded operator e~~^^, and the vector 

|A) := ^^^ e n (28) 

||e-^^MO)|| 

defines a (^, a^)-KMS state on A — Aiy ® Ay^r and it satisfies 

Lx\X) = 0. (29) 

Proof. To show that |0) G Doni(e^"-^^) we check that the Dyson series 

^(-A)« r\t,--- r""di„a[',J/)...<(/)|0) (30) 

n>o -^0 -^0 
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converges. We write the interaction operator conveniently as / = GQ[g], where 
G = A + A^ and g{x) = aS{x'^)x^p{x^,{x)) has support in Wr, c.f. itlfill . In l|30ll 
we have set, for real s, 

To see that e~^^^Q[g]e''^^ is well defined for < s < n/a one shows that since 
g is supported in the right wedge, the map t i-^ &'^^^ Q[g]e~^^^^ = Q[g o B_t] 
has an analytic continuation into the strip < Imi < n/a, and it is continuous 
at the boundary of the strip (i S M, t S i-R). This argument is actually part of 
the proof of the Bisognano-Wichmann theorem, [BiWi|. It follows in particular 
that the integrals in 1(3(1)1 are well defined and that furthermore 



sup 

0<lm s<7T /a 



.Ul){N + l)-'/' 



= C < oo, 



where N is the number operator on Fock space. Since |0) is the vacuum on the 
field part, and each interaction term aP,(/) can increase the particle number by 
at most one we have the bound ||a?(^ (/)••• a?j^(/)|0)|| <C"Vn^.. It follows that 

the series (13011 converges (for all values of A) and hence |0) G Dom(e^'^^^). 

The facts that |A) defines a (^,a^)-KMS state and that Lx\0) = follow 
from Araki's perturbation theory of KMS states, and from perturbation theory 



of standard Liouville operators, see |D.TP| . ■ 

We are now in a position to state the full result, of which Theorem 12. 21 is an 
immediate consequence: 

Theorem 3.5 Assume that the Fermi Golden Rule Condition ^^ is satisfied. 
There exists Xq so that for all < |A| < Aq, for all density matrices g onH and 
for all B e A 

lim Tr ga^iB) = {\\B\X). 

a — >oo 

Proof. We show in Section 13.31 that the result follows if the spectrum of L\ 
is purely absolutely continuous with the exception of a single simple eigenvalue 
at zero. These spectral characteristics are shown in Theorem 14.21 ■ 

3.3 Reduction to a spectral problem 

We reduce the proof of Theorem 13.51 to a spectral problem via the following 
simple lemma, which is a variant of the Riemann-Lebesgue lemma: 

Lemma 3.6 Let Ti he a Hilbert space, 4> <E Ti, A a subalgebra of B{T-L) whose 
commutant we denote by A' , and let L be a self-adjoint operator on Ti. Suppose 
that A'(j> is dense in Ti, that e^^'^ Ae~^^'^ C A, Vt, that Lip = 0, and that on the 
orthogonal complement of 4>, L has purely absolutely continuous spectrum. 
Then we have 

lim Tr g e'^^Be"'^^ = (0,^0), (31) 

r — >oo 

for all A £ A and for all density matrices Q < g <E C^{T-L), Trg = 1. 
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Proof. We may diagonalize g = X]^i-PnlV'n)(V'n|, where ipn & Ti- and the 
probabiHties < p„ < 1 sum up to one. So it suffices to show (I31|l for a 
rank-one density matrix g = |V')(^I- Given any e > there is a B' G A' s.t. 
Wijj — B'(j)\\ < e. Thus by replacing i/; by B'<j), commuting B' and e'-^'^^e^'^'^, 
and by using the invariance of (j) under e"'^"^ we obtain 

Trge'^^Be^'^^ = (V', e'^"Be-'^"V) = {^, B'e'^^B^) + 0(e), (32) 

where the remainder is estimated uniformly in r. Since the spectrum of L is 
absolutely continuous except for a simple eigenvalue at zero with eigenvector 
(j), the propagator e'^"^ converges in the weak sense to the rank-one projection 
\(j}){(j)\, as r ^ oo. Using this in ll32ll . together with the facts that {ij^tB'^i) = 
1 -|- 0(e), and that e can be chosen arbitrarily small yields relation ()31|l . ■ 

We apply Lemma l3.(il with L = L\ and 4> = \^)- The density of A!\\) 
follows from the KMS property of |A), the invariance of A under e'''^^'^ • e"'^^"^ 
follows from Lemma l331 and the relation La|A) = is shown in Proposition |OI 
It remains to prove that on the orthogonal complement of |A), L\ has purely 
absolutely continuous spectrum. 

4 Spectral analysis of L\ 

The spectrum of the operator Ld consists of two simple eigenvalues ±E (eigen- 
vectors |±,=f)) and a doubly degenerate eigenvalue at (eigenvectors |±,±)). 
Lp has absolutely continuous spectrum covering the entire real axis, and a single 
embedded eigenvalue at the origin. This eigenvalue is simple and has eigenvector 
|0,F). It follows that Lq has absolutely continuous spectrum covering the axis 
and three embedded eigenvalues at 0, ±E, the one at being doubly degenerate. 
Our goal is to show that the nonzero eigenvalues are unstable under the 
perturbation A(/ — JIJ), and that the degeneracy of the eigenvalue zero is 
lifted. We do this via spectral deformation theory, showing that the unstable 
(parts of the) eigenvalues turn into resonances located in the lower complex 
plane. 

4.1 Spectral deformation 

For the spectral analysis it is useful to consider the unitarily transformed Hilbert 
space L'^{W^, d>cd'^~^kj_) of one-particle wave functions of the field, determined 
by L^{R'^,d''x) 3 f^Wf with 



iWf){M:,k^) := \/lu± cosh }€ f^Lu^smhti^kj^), (33) 



where uj± := \/\E±\^ + rn^ and where / is the Fourier transform of /. The 
advantage of this representation of the Hilbert space is that the operator K, 
defined in 10, takes the particularly simple form K = id^. The transformation 
W lifts to Fock space in the usual way. We do not introduce new names for 
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spaces and operators in the transformed system. The Liouville operator Ijl6|l is 

La = Ld+Lf + \V, 

Lq = Lu + aLp, Ld = ffD«)]l2-]l2«)LfD, LF = dT{id^), (34) 

V = I~JIJ, I = G®t2®^{a\g)+a{g)) (35) 

v2 

acting on the Hilbert space Ti — 'C'^®C?®!F, where J- is the bosonic Fock space 
over L^(IR'', dx: d'^~^fcj^). In (1^ G is the 2x2 matrix with on the diagonal and 
1 on the off-diagonals, and g{}<,k^) = {W^~^/'^xip{x^\r=o)){>c,k^_) is given in 
l|n||. The action of jp, (EHJ, is given by {jff){>c,k^) = f{-K,k^_). 
We describe now the complex deformation. Let e M. The map 

Mki, ...,Kn):= {Ue^) [ku-.., ^n) ■■= e"'(-i+--")V(>^i, • • ■ , ^n) 

defines a unitary group on T (we are not displaying the variables kj_ in the 
argument of tp since Ug does not act on them) . An easy calculation shows that 

Lxie):=UgLxU;=Lo{e) + XV{0), Lo{9)=Lo-a0N, V{9) = 1(6) - JI{e)J, 

(36) 
where A^ = dr(]l) is the number operator on T and 

1(9) = G®h®-^ {a^(e''"ff) + a(e'^"5)} , (37) 

where we have put the complex conjugate 9 in the argument of the annihilation 
operator in 113 7|l in view of the complexification of 9. 

Lemma 4.1 Let 

9o{m,d):^l Z t'^^y'ttzl (38) 

where m > is the mass of the field and d is the spatial dimension. We have 
e'^'^Wn-^^^h € L'^iW^.dKdk^;) for all 9 e C satisfying \9\ < 9q and for all 
h £ S{M.'^,C). Moreover, for \9\ < 9o, L\{9) is a closed operator on the dense 
domain V = Dom(Lo) n Dom(Af). 

Proof. Lo{9) is a normal operator, so it is closed. Assume we know that 
e'^^VFO-i/25(R'',C) C L"^ {W^ , dxdk^} , and recall that a;i/9(x,|r=o) G SiW^.C)- 
Then, for ImS ^ the perturbation V{9) is infinitesimally small w.r.t. Lq{9), 
so Lx{9) is closed by stability of closedness. For Im0 — the operator L\{9) is 
even selfadjoint. 

Let h£S{W^X). According to ^ we have 

(W^rj-^/2/i)(>f,fc^) =/J(cJiSinh>f,fc^). 
Since h £ S we have that for any integer n there is a constant C„ s.t. 

h{uj± sinh>f, fcj 



Cn 
< 



1 + [m^ sinh^ >c + [fc^p cosh^ . 
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For 7TT- = we thus obtain (using an obvious change of variables) the estimate 

r I. . f - 2 _ /• p2e'|3<l 
/d^e^^l^l/ dfc, h(cu±smh>c,k.) < Cn d>c — ^— p (39) 

Jm V-1 Jr [cosh>*]'i-i 

which is finite provided 9' = \lme\ < {d - l)/2. If m 7^ then the l.h.s. of ^ 
is bounded from above by 

/nt2 



[l/2 + m2sinh^>f]"[l/2+ |fc_L|2]» 
which is finite if 0' = \Im9\ < n, and n can be chosen arbitrarily large. ■ 

4.2 Spectra of Lx{0) and of La 

The goal of this section is to prove the following result. 

Theorem 4.2 Suppose the Fermi Golden Rule Condition Ufi\} holds. There 
is a Xq > s.t. «/ < |A| < Aq then the spectrum of L\ consists of a simple 
eigenvalue at zero and is purely absolutely continuous on the real axis otherwise. 

Remark that the spectrum of Lq{9) — L^ — aON consists of the isolated 
eigenvalues ±E (simple) and (doubly degenerate), and of the fines of con- 
tinuous spectrum {M — ianlm9}n=i,2,...- We now analyze the behaviour of the 
eigenvalues of Lo{9) under the perturbation XV (9). The strategy is to show that 
the eigenvalues zLE are unstable under the perturbation, and that the degener- 
acy of the eigenvalue zero is lifted. Note that the kernel of L\ is non-empty by 
construction, see ll29|l . 

Proof of Theorem \4-S\ The central part of the proof is the control of 
the resonances bifurcating out of the eigenvalues zLE and 0, see Lemma 14.31 
A standard analyticity argument then implies Theorem 14.21 The latter can 
be summarized as follows: one checks that for all complex z with Imz > 0, 
9 I— > {tj^,{L\{9) — z)~^4>0) is analytic in < |6'| < ^o, Im^ > 0, and continu- 
ous as \uy9 [ 0, for a dense set of deformation analytic vectors 1/', 4> (take e.g. 
finite-particle vectors of Fock space built from test functions /(><) with com- 
pact support). As is well known, the real eigenvalues of L\{9) coincide with 
those of L\ , and away from eigenvalues the spectrum of L\ is purely absolutely 
continuous. 

We now present in more detail the resonance theory. 

Lemma 4.3 Let 9' = Im9 > 0. There is a Ai (independent of 9) s.t. if 
|A| < Aimin(l,6'') then in the half-plane {Imz > —9' /2\ the spectrum of L\(9) 
consists of four eigenvalues £±(A) and £q{X) and only (which do not depend 
on 9). 

Moerover, we have ej{X) — Cj — X'^e, + o(A^), where j = +, — , 0, and 

4'^ =iIm4''-i(l + e"'"»)C, (40) 

whith £, given in 
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We remark that the Fermi Golden Rule Condition IjlSII asserts that ^ > 0. 

Proof of Lemma \4-3\ By an argument of stability of the spectrum it is not 
difficult to show that the spectrum in the indicated half plane consists of four 
eigenvalues only. 

The position (at second order in A) is governed by so-called level shift oper- 
ators, see e.g. |M2| and references therein. We explain this with the help of the 
Feshbach map |BaFSj . 

Let e be an eigenvalue of Lo{9) and denote the corresponding (orthogonal) 
eigenprojection by Qe = PePvac, where Pe is the spectral projection of L-q onto 
e and Pvac projects onto the vacuum in !F . Set Qg := 1 — Qe and denote by 
X — Q^XQ^ IkhiiQ ^^^ restriction of an operator X to RanQ. A standard 
estimate using Neumann series shows the following fact. 

Lemma 4.4 There is a constant A2 (independent of 6) s.t. if |A| < A2 min(i<^, 9') 
then, for each eigenvalue e of Lq, the open ball of radius 9' /2 arounde, B{e, 9' /2), 
belongs to the resolvent set of L^{9). 

It follows from Lemma ^3| that the Feshbach map 

Fe^ALx(O)) := Qe (e - \^V {9)Q ,{V^{9) - z)-^Q,V{9)) Qe (41) 

is well defined for all z S B{e,9' /2). This map has the following remarkable 
isospectrality property |BaFS| : for all z e B{e,9' /2), 

z £ spec{Lx{9)) ^^ze apec{FeALx{9))) . (42) 

Thus it suffices to examine the spectrum of the operator Fe^z{L\{9)) which acts 
on the finite dimensional space RanQe ■ We expand the resolvent in l|lTll around 
A = and consider spectral parameters z ^ e + 0(A) to obtain 

FeALx{9)) = Qe{e- X^V{9)Q^(Lo{9) - e)-^Q^V{9)) Qe + oiX""), 

where hmA^o o(A^)/A^ = 0. We now use analyticity in 9 to conclude that 

Fe,,{Lxm = Qe{e- X^VQeiLo - e - iO+)-'QeV) Qe + o(A2), (43) 

where iO+ stands for the limit of ie as e | 0. The operators 

Ae := QeVQeiLo ~ 6 - iO+)-'QeVQe 

are called level shift operators. For e = ±E they reduce in the present case 
simply to numbers (dimRanQe — 1), while Aq corresponds here to a 2 x 2 
matrix. Using the expression IJ35II for V one can calculate explicitly the level 
shift operators (see also |BaFSI IMll IM2| for more detail on expHcit calculations 
in related models). 

Lemma 4.5 In the basis {| — ,— ), |+,+)} o/RanQo we have 



Ao = i^e-"- 
where ^ is given in \21]) 



&-"- -1 
-1 e'^f 



and ImA±B = (l + e ^^)i 
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Remark. The Gibbs state of the detector at inverse temperature /3 = 27r/a 
(represented by a vector ex [l,e~'^^'°]) spans the kernel of Aq. 

This lemma together with lll3|l and the isospectrality l(l2ll shows the expan- 
sions ijifijl and IpU- This proves Lemma ^31 and at the same token, concludes 
the proof of Theorem 14.21 ■ 

5 Proofs of Proposition 12.11 and of Lemma 13.31 

Proof of Provosition \2.1[ The coupled Liouville opertor Ill6|l has the form L\ = 
Lq + XI, where / = GQ[g\ with G = A + A'' and g{x) = aS{x'^)x^p{x^,(x)) has 
support in Wr. Essential selfadjointness of L\ can easily be shown using the 
Glimm-Jaffe-Nelson commutator theorem, see e.g. [FMj, Section 3. 
The Araki-Dyson series expansion gives (weakly on a dense set) 

e"^-Me-"^- = ^A" / dti / di2--- / dtjG{h)Q[g o B^at,], ■■■ 

„-n "'0 Jti Jt„-i L L 

■]], (44) 

where we set G{t) = e'*^°Ge"'*-^°. For M £ A any element M' e A' com- 
mutes termwise with the series (011, hence M'e^*^^Me~^'^^^ = e'*''"^Me"'*^^M'. 
Therefore we have e'*''^^^e"'*^^ & A. ■ 

Proof of Lemma \S.Sl Essential selfadjointness is shown using the Glimm- 
Jaffe-Nelson commutator theorem, see e.g. [FMJ Section 3.1. The fact that L\ 
and L\ define the same dynamics on A is easily derived by using that L\ — L\ 
belongs to the commutant of ^ (and e.g. applying the Trotter product formula), 
see also IFMI. ■ 



[G(t„)g[goB_atJ,e'*^«Me-"^° 
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